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Abstract. - We propose an experimentally relevant scheme to create stable solitons in a 
three-dimensional Bose-Einstein condensate confined by a one-dimensional optical lattice, us- 
ing temporal modulation of the scattering length (through ac magnetic field tuned close to the 
Feshbach resonance). Another physical interpretation is a possibility to create stable 3D "light 
bullets" in an optical medium with a longitudinal alternating self-focusing/defocusing struc- 
ture, and periodic modulation of the refractive index in a transverse direction. We develop a 
variational approximation to identify a stability region in the parametric space, and verify the 
existence of stable breathing solitons in direct simulations. Both methods reveal that stable 
solitons may be supported if the average value of the nonlinear coefficient (whose sign corre- 
sponds to attraction between atoms) and the lattice's strength exceed well-defined minimum 
values. Stable localized patterns may feature a multi-cell structure. 



The creation of the Bose-Einstein condensates (BEC) in vapors of alkali metals has opened 
an opportunity to investigate nonlinear interactions of atomic matter waves. One of the most 
fundamental aspects of these studies is generation of solitons. Dark solitons had been created 
in effectively one-dimensional (ID) BEC with repulsive interatomic collisions [1]. Then, bright 
solitons were observed in ID BECs with attractive collisions [2,3]. An issue of obvious interest 
is to develop methods for control of the matter-wave solitons. A promising approach consists 
in varying the scattering length (SL) of the interatomic collisions by means of an external 
magnetic field through the zero-SL point close to the Feshbach resonance [4]. In particular, 
it was shown that an abrupt change of the SL can result in splitting of a soliton into a set of 
secondary ones [5]. 

The application of ac magnetic field may induce a periodic modulation of the SL, opening a 
way to the "Feshbach-resonance management" (FRM) [6] . Quite a noteworthy FRM- induced 
effect is a possibility to create self-trapped oscillating BEC solitons (breathers) without any 
external trap. The underlying mechanism is balance between compression and expansion 
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cycles, corresponding to the time intervals when the SL is, respectively, negative and positive. 
In fact, this possibility was first realized in terms of the transverse light propagation in a 
bulk layered nonlinear medium, where the Kerr coefficient alternates between positive and 
negative values [7]. The respective model is formally equivalent to the two-dimensional (2D) 
Gross-Pitaevskii equation (GPE), see Eq. (JJJ below, in which the SL is periodically modulated 
in "time" (analog of the propagation distance in [7]) not harmonically, but as a piecewise- 
constant function. Then, the BEC model based on the GPE with the harmonic modulation 
of the SL was directly investigated in Refs. [8-10], with a conclusion that the FRM makes 
it possible to stabilize, without the use of an external trap, 2D breathers, but not the 3D 
ones. The absence of the stabilization in the 3D case was also concluded in the optical model 
considered in Ref. [7], where the 3D soliton would correspond to a spatiotemporal "light 
bullet". 

A remarkable feature of these results is that the 2D breather is stable despite a possibility 
of the collapse in the 2D GPE with a negative SL. On the other hand, it has been recently 
demonstrated that solitons in BECs with a constant negative SL may be completely stabilized 
by means of an optical lattice (OL), i.e., a spatially periodic potential. It is true in both 2D 
and 3D cases [11]. Moreover, it was also shown that the 2D and 3D BEC solitons can be 
stabilized by low- dimensional OLs, i.e., respectively, by a quasi-lD (Q1D) lattice in the 2D 
case, and by a quasi-2D (Q2D) lattice in the 3D space [12]. However, 3D solitons cannot be 
made stable in a Q1D OL. 

These results suggest a question whether 3D solitons can be stabilized by a combination 
of a Q1D lattice and FRM. The issue has practical relevance, as a Q1D OL can be easily 
created, illuminating the BEC by a pair of counterpropagating laser beams that form a periodic 
interference pattern [13]. This problem is different from that considered in Ref. [10], where 
stabilization of a localized 3D structure was provided by adding a tight ID parabolic trap to 
the model, which made it nearly two-dimensional. 

In this paper, we demonstrate that the combined OL-FRM stabilization of 3D solitons is 
possible indeed. First, we introduce the model and propose an experimental scheme. Then, 
we predict stability conditions for the 3D solitons by means of a variational approximation 
(VA). Finally, we demonstrate the soliton's stability in direct simulations. We also demon- 
strate that these results find another physical realization in nonlinear optics: a possibility to 
create a stable 3D "light bullet" (spatiotemporal soliton) in a bulk medium which combines 
an alternating self-focusing/defocusing structure in the longitudinal direction, and periodic 
modulation of the refractive index in one transverse direction. 

The model. The GPE for the single-particle wave function ip, including a time- 
dependent (FRM-controlled) nonlinear coefficient g(t) and an external potential V(r,t), in 
normalized units is 

nPt= [-(1/2)V 2 + V(r,t) + g(tM 2 ] 4, . (1) 

We start by considering a BEC in the ground state of a radial (2D) parabolic trap with the 
(time-dependent) frequency u>±(t), supplemented, in the longitudinal direction, by "end caps" , 
induced by transverse light sheets. The configuration is like the one used to create soliton 
trains in the Li 7 condensate [3]. Then, a ID lattice in the axial direction, whose period is 
normalized to be tt, is adiabatically switched on, by increasing its strength e from at t = 
to a final value Sf at t = see fig. ^ Thus, the full potential is 

V(t, t) = e(t) [1 - cos(2z)] + (\/2)ul{t) B 2 + V (z) , (2) 

where g is the radial variable in the plane transverse to z, and the axial end-cap potential 
Vq(z) is approximated by a deep and wide potential box. 
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Fig. 1 - The time dependence of the nonlinear coefficient, g, transverse-confining frequency, ui±, and 
optical-lattice strength, e, in the numerical experiment which leads to the establishment of stable 
3D breathing solitons supported by the combination of the quasi- ID lattice and Feshbach- resonance 
management (FRM). The shaded area indicates rapid oscillations of g, which account for the FRM. 



At some moment t = t\ < £2, we begin to linearly decrease an initially positive nonlinear 
coefficient g(t). It vanishes at the moment t = t%, and remains zero up to t = £3, when we 
start to gradually switch on the rapid FRM modulation of the nonlinear coefficient (in fig. ^ 
g(t) is denoted as go(t), up to £ = £3). At t > £ 4 , <?(£) oscillates with a constant amplitude gn 
around a negative average value #of > 

g(t) = 9of + 9if sin(fi£) . (3) 

Finally, the radial confinement is switched off within the time interval £3 < £ < £4. After that, 
a soliton, if any, may only be supported by the combination of the Q1D lattice and FRM 
(the axial width of the established soliton is assumed to be much smaller than the distance 
between the end caps, hence they do not affect it). 

Numerical experiments following the path outlined in hg. ^ make it possible to create 
stable 3D solitons. Before showing the results, we first resort to the VA, in order to predict 
conditions on the strength of the OL and size of the negative average nonlinear coefficient got, 
which are necessary to support 3D solitons. 

Note that the eventual form of Eq. QJ, i.e., with g and V taken, respectively, as in Eqs. (J3J) 
and J5J, where e is constant and u>± — Vo(z) — 0, has an alternative physical interpretation in 
terms of nonlinear optics. Upon interchanging variables £ and coordinates z (chosen to be a 
propagation direction) Eq. (JJJ describes pulse propagation in a bulk (3D) medium composed of 
alternating self- focusing and self-defocusing layers, similar to that introduced in Ref. [7], with 
an additional periodic modulation of the refractive index in the transverse direction x. Stable 
3D solitons in the present model imply the existence of stable fully three-dimensional spa- 
tiotemporal solitons ("light bullets") in this medium. Thus far, only quasi-2D spatiotemporal 
solitons were created experimentally in bulk media [14]. 

Variational approximation. The VA was successfully applied to the description of 
BEC dynamics under diverse circumstances, see, e.g., Refs. [5,8-12,15] and a review [16]. 
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Equation is derived from the Lagrangian 

L = w [ [ityW-W) ~ \<l> e \ 2 ~ \<Pz\ 2 ~ 9(t)M 4 ~ 2V\i } \ 2 ]gdgdz. (4) 



As it is frequently done, we choose a variational ansatz for the solution based on the complex 
Gaussian with an amplitude A(t), radial and axial widths W(t) and V(t) respectively, b(t) 
and j3(t) being the corresponding "chirps" 

i/)(r,t) = Aexp(-Q 2 [l/(2W 2 )+ib] -z 2 [l/(2V 2 )+ip]) , (5) 

The reduced Lagrangian is obtained by inserting the ansatz (JSJ) into Eq. (JIJ. The equations 
obtained by varying the reduced Lagrangian yield the conservation of the number of atoms, 
E = 7r~ 3 / 2 J |^| 2 dr = A 2 W 2 V , and dynamical equations for the widths, 

W = -L-u 2 ± (t)W+ f. 9 ® , (6) 

V = ± - AeV exp (-V 2 ) + ^ . (7) 

A necessary condition for the existence of a 3D soliton in the present model can be derived 
from these equations in a crude approximation, cf. Ref. [8]. To this end, we assume that 
gn in Eq. © is small, while f2 is large. It is also conjectured that the average value W of 
the soliton's radial size, W, is large (see below). Further, in the lowest approximation, the 
soliton's size in the axial direction may be assumed constant, V(t) w Vb, as determined by 
the relation 

4 £ V/ 4 exp(-V/ 2 ) = 1, (8) 

that follows from Eq. (JSJ) where the last small term (~ W~ 2 ) is dropped. Equation ijSJl has 
real solutions if the OL strength exceeds a minimum (threshold) value, 

£>ethr = e 2 /16«0.46. (9) 

At e > £thrj Eq- © has two real solutions, which implies the existence of two different solitons. 
It seems very plausible (cf. Refs. [11,12]) that the narrower one, corresponding to smaller Vo, 
is stable, and the other one is unstable. 

Next, replacing V(t) by Vb in Eq. JJJ) (we set cu± — 0, to consider the possibility of the 
existence of the 3D soliton without the radial confinement) and substituting git) from Eq. 
we look for a solution as W(t) w W+Wi sin (fit). For large O, the variable part of the equation 

yields W\ = —Egn/ (^/8W 3 Vo^l 2 ^j ■ Then, the consideration of the constant part of Eq. |JBJ|, 
with regard to the first correction generated by the product of the oscillating terms in W(t) 
and g(t), produces a result, 



w 4 3 (Em 



\ 2 -1 

J (JS IflQf] - V8V&) . (10) 



4-s/2Vb V ^ 

An essential corollary of Eq. (|10fl is that a necessary condition for the existence of the 3D 
soliton is (recall g^i is negative) 

Iffofl >(|<?of|) min = V8V /£. (11) 



Direct simulations presented in the next section show that this condition holds indeed, albeit 
approximately. Note that the above approximation, which assumed large W , is valid only 
when |gof| slightly exceeds the minimum value defined in Eq. (|1 lfl . 
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Fig. 3 



Fig. 2 - Evolution of the absolute value of the wave function in the pattern comprising five cells of the 
optical lattice. The normalized parameters are <?of = 22.5, gu = 4<? f, £f = 25, Q = 40, w±(0) = 0.3, 
ti = 30, t a = 100, t 3 = 120, and t 4 = 130. Snapshots are taken at t = (a), 20 (b), 50 (c) 110 (d), 
200 (e), 2000 (f). 

Fig. 3 - Evolution of the central peak amplitude from fig. [5] 



Numerical results. - We simulated both the full GPE (Eq.Q using an axisymmetric code, 
and the variational equations and |7|> ■ In both cases we used Runge-Kutta methods. In the 
former case, the simulations followed the path outlined in fig. ^ Example of the numerical 
results, which is quite generic for strong lattice, is displayed in figs. [5] and [21 For instance, for 
87 Rb atoms, the OL period is A = 1.5 /im, the FRM and initial radial-confinement frequencies 
are, respectively, = 132 kHz and w_l(0) = 990 Hz, the lattice depth is e = 25 recoil energies, 
and the effective nonlinear coefficient is na = ± 4 • 10~ 6 m, where n is the number of atoms 
per the lattice cell, the total number of atoms being in the range of 10 4 — 10 6 . The respective 
values of the normalized parameters are given in fig. [3 caption. 

In fig. [3 we show 3D snapshots of the wave-function pattern occupying five lattice cells, 
taken at times corresponding to various stages of the numerical experiment, cf. fig-H Compar- 
ison of the last two snapshots, (e) and (f), shows that condensate develops a robust structure, 
which remains unchanged over many FRM cycles. 

Figure [3] displays evolution of the central-peak's amplitude in the same pattern. After an 
initial transient, the stable structure is established, featuring breathings without any system- 
atic decay. 

Actually, the stable pattern like the one shown in fig. [21 is a set of uncoupled fundamental 
solitons, each being trapped in a single lattice cell. This conclusion follows from additional 
numerical experiments, in which removal of the atomic population from any subset of the cells 
did not, in any tangible way, affect the localized states in other cells; in particular, the soliton 
sitting in a single cell is as stable as any multi-cell pattern. Thus, the conclusion is that the 
FRM may support a stable 3D soliton (breather) confined to a single cell of the Q1D lattice, 
and the solitons trapped in adjacent cells do not interact. In fact, this conclusion complies 
with the variational ansatz l|5"jl. which ignores any possible spatial oscillations induced by 
the OL, hence it indeed implies a soliton essentially confined to a single cell (cf. the VA for 
stationary models with the multi-dimensional OLs developed in Refs. [11] and [12]). 

We have collected results of systematic GPE simulations and compared them with pre- 
dictions of the VA based on simulations of Eqs. (HO and 0) , as shown in fig. 01 As is seen in 
fig-GI a )j on ly the bottom part (corresponding to smaller values of the FRM frequency f2) of 
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Fig. 4 - Stability regions for the 3D solitons as predicted by the variational approximation, and 
found from direct simulations of the Gross-Pitaevskii equation in (a) the (got i fi) plane, and (b) the 
(e, fi) plane. Other parameters are as in fig. [5] except that <?of = 30 in (b). The vertical line in (a) 
corresponds to the minimum value of \gof\ predicted by Eq. IIH . 



the VA-predicted stability region in the (gof, Q) plane actually supports stable 3D solitons. 
The variational estimate Hll|) for the minimum size of the average nonlinear coefficient nec- 
essary for the existence of the 3D soliton in the Q1D lattice is borne out by the simulations, 
although approximately. 

Lastly, the simulations definitely demonstrate the existence of a minimum of the OL 
strength e which is necessary to support the 3D solitons, as it was predicted by the VA, 
see Eq. (0. Figure Efb) shows regions of stability in the (e, f2) plane. From this figure we 
conclude that VA predicts stability regions quite well, but the numerical value of the e t hr 
is somehow larger than that predicted by the simple approximation given by Eq. For 
simulations presented in figs. |2Jand|21we have chosen the value of e well above the threshold. 

Conclusions. In this work, we have proposed a scheme to stabilize 3D solitons in 
BECs where the interatomic interaction is attractive, on average. The scheme is based on 
the combination of the Feshbach-resonance management and quasi-lD optical lattice. In 
previous works, it has been shown that these two methods applied separately can stabilize 2D 
solitons, but not 3D ones. The combined method proposed here can be easily implemented 
in the experiment, and thus opens a way to the creation of 3D solitons in BECs. Another 
physical implication of the results is a stable 3D "light bullet" in an optical medium with 
an alternating longitudinal self-focusing/defocusing structure, and periodic modulation of the 
refractive index in a transverse direction. 

The possibility to achieve the result was demonstrated within the framework of the vari- 
ational approximation (VA) and verified in direct simulations. In particular, the VA predicts 
that the OL's strength and average value of the oscillating nonlinear coefficient must ex- 
ceed certain minimum values. These predictions are indeed corroborated by the simulations. 
However, the full stability region predicted by the VA extends much farther into the high- 
frequency region than the actual stability area found in the direct simulations. The stable 
patterns found in the simulations may feature a multi-cell structure, which in the case studied 
here form a set of non- interacting fundamental solitons, each trapped in a single cell of the 
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OL, in good qualitative agreement with the description provided by the VA. In weak lattices 
one may expect to observe interaction between individual solitons. This will be a subject of 
further investigation. 
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